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Abstract 

We discuss non perturbative corrections to the Coulomb-like potential of heavy 
quarks at short distances. We consider both the standard framework provided by 
infrared renormalons and the assumption that confinement does not allow weak fields 
to penetrate the vacuum. In the former case the leading correction at short distances 
turns out to be quadratic in r for static quarks. In the latter case we find a potential 
which is proportional to r as r — > 0. We point out that similar effects arise due to a 
new kind of non perturbative correction proportional to 1/Q 2 , which is unaccounted for 
by the operator product expansion and which was recently discussed within a different 
framework. Phenomenological implications of the linear correction to the potential are 
briefly reviewed. 
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1. Analytic studies of non-perturbative effects in QCD are limited by the fact 
that there does not exist a systematic procedure to account for confinement. We may 
attempt, however, to sort out and to parametrize its effects by considering the case 
when the perturbative part is dominant and the infrared sensitive contributions appear 
only as small power corrections in a large mass scale Q. Well known examples of this 
kind are provided by the applications of the operator product expansion such as the 
QCD sum rules More recently, renormalons have emerged as a more universal 
tool for enumerating the power-like corrections even in situations where the operator 
product expansion does not apply ( for furthur references and review, see e.g., H). 

In this note we will apply these techniques to the short- distance interaction between 
heavy quarks, Q,Q. As is well known, this interaction is dominated by a Coulomb-like 
potential: 

IimV(r) =-° Fa ^ (1) 
r^O r 

where the coefficient Cp = 4/3 refers to the color-singlet channel. We will argue 
below that the character of the leading power corrections depend on the properties of 
confinement. In particular, the standard renormalon technique implies the effect of 
confinement to become important only at distances R cr ~ ^qcd an d the leading power 
correction to (P turns to be: 

\im5V{r) renormalon = c 2 r 2 . (2) 

1 >0 

Similar results for the static potential have already been obtained in the literature by 
various techniques @, ||, |(| . 

Although infrared renormalons probably correctly indicate the presence of some of 
the power corrections, it is an open questions whether they exhaust all the possibilities. 
In case of the static potential, infrared renormalons do not give any indication to a 
formation of a string at large distances and could miss, therefore, the leading correction 
at short distances as well. In fact, we will consider the power corrections to the static 
potential (|1|) in an alternative model which implements the assumption that the effects 
of confinement become important once the color fields become weak. In such a case 
(see below ) we find in contrast to Eq. (^), a leading correction which is linear in r: 

lim SV (^non-standard = C\T (3) 

A crucial element to ensure such corrections appears to be the existence of small size 
non perturbative fluctuations in QCD. 

Power corrections which go beyond the OPE and infrared renormalons have recently 
been highlighted both within a general dispersion approach |7|, || and the ultraviolet- 
renormalon technique f9|. Most remarkably, indications from lattice simulations for the 
existence of such corrections have very recently been obtained in ||10[| . More specifically, 
one argues that the effective coupling may contain terms proportional to 1/Q 2 '- 

*eff(Q 2 ) = a^XQ 2 ) + c Q2 ^ (4) 

while infrared renormalons only indicate the presence of corrections proportional to 
1/Q A , corresponding to the gluon condensate. 



1 



It is amusing to observe that the 1/Q 2 piece in (Q) also would give rise to a linear 
correction to the potential. Thus it is interesting to speculate that the mechanism of 
the repulsion of weak fields by the vacuum could be a possible mechanism responsible 
for the correction @ if its existence [plj is confirmed. We would like to emphasize 
that although the static potential is not directly observable, it can be measured on the 
lattice (for a recent study including further references see |l ]]]). In particular, although 
the existence of a linear potential at large distances is well established and there is no 
indication to its change to the r 2 behaviour as predicted by (|2|), much more precise data 
is needed to either confirm or reject the presence of a linear term at short distances. 
On the other hand the linear correction Eq. (|3]) would also affect the by-now-standard 
theory of the bound states of heavy quarks (for a review and further references see, e.g., 
[|12j) and we will address this issue as well. 

2. To evaluate the renormalon contribution to V(r) consider first the one gluon 
exchange potential 

V(r) = -C F |ci 3 k(4vra s (k 2 )) eXp( g- r) (5) 

where the leading logarithmic corrections have been incorporated into the running cou- 
pling constant, a s (k 2 ), and the fact that for static quarks it depends on k 2 is crucial 
for further argument. Beyond the leading logs, Eq. (|5|) holds in an Abelian case and in 
the limit of large Nf in QCD ||. To find the renormalon contribution we write 

whee b is the first coefficient in the expansion of the /3-function (for simplicity we 
confine ourselves to a one-loop /^-function), and Aqcd is the position of the Landau 
pole. 

Next, we substitute (^) into (Q) and perform the integration over directions of k to 

get 

V(r)= -%r*r4£M( * V"* ,7) 
3 Jo Jo kr \Aqcd J 

Moreover, and since we are interested in V(r) at small r, we expand in kr: 

V (r) = - I|l J™ Affi D da J dk{k- 2 ^ + k 2 - 2 ^r 2 + ...), (8) 

The renormalon poles at a = l/2&o, 3/2&0-- are clearly seen in the above and this 
suggests a modification of the Coulomb potential at small r of the form: 

6V(r) = coAqcd + c 2 A 3 QCD r 2 + ... (9) 

where Co t 2 are constants and it is noteworthy that they do not contain any small cou- 
pling. This has already been noted in 0. Since we expand in kr and k ~ Aqcd at 
the renormalon poles, Eq. (H) is valid at distances small compared to the confinement 
radius: 

r « Aqcd- (10) 
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This result has been obtained via a Borel transformation to conform to the current 
usages in renormalon analysis. It may also be seen by noting that the pole at k — Aqcd 
in the standard, leading log expression for a s , 

6 log k 2 /A z QCD 

implies that the Fourier integral above is undefined. According to the general theory 
of singular integrals, we may write 

a - (k2) - pp + icA « CD<S(fc2 " A «™»' (12> 

where c is an arbitrary constant. When integrated the last term, ick^ CD 5{k 2 — Aq CD ), 
is easily seen to give a correction to the potential, 

5V(r) = c Sin(rAQCT) (13) 
r 

which, upon expansion, reproduces the type of corrections given in Eq. ([|). 

Eq. (||) remains qualitatively unchanged upon inclusion of higher loop corrections 
as well. The results of explicit calculations (which are known now up to three loops 
fll3"| ) are usually represented as an expansion, in ajjg'. 

V(k 2 ) = -C F - \ 2 MsK . (14) 

It is straightforward to see that renormalon contribution associated with higher powers 
of a(k 2 ) give the same type of power corrections as above. Indeed, in the approximation 
of the one-loop /3-function, 

a 2 (k 2 ) = ^A QCD -^-a(k 2 ). (15) 
2b A-QCD 

Applying the differentiation with respect to Aq C d at the last step, i.e. to Eqs. fl9],|13"D we 
immediately see that the Fourier transform of the a 2 (k 2 ) term brings the same kind of 
power corrections as above. This is true for higher powers of a(k 2 ) as well so that Eq. (|9|) 
is reproduced by renormalons associated with any order in a(k 2 ) in Eq. fll4|). Finally, 
inclusion of higher loops in the /^-function itself is known [14]] to modify the renormalon 
poles to renormalon cuts, the position of the singularity remaining unchanged. This 
would bring powers of logs in the correction to the potential but would not change the 
power of r, r 2 . 

It is worth pointing out that the renormalon contribution could be tried also at large 
distances. The resulting potential does not reproduce the linear rise at large r and is 
readily seen to be devoid of physical meaning. The reason is that renormalons are a pure 
perturbative construct and do not produce any hint that at large r a string is formed. 
This failure demonstrates a limitation of the renormalon technique - the existence of 
the Landau pole in the infrared region does not automatically imply confinement. 

To go beyond the perturbative expansion ([14] ) for V(r) in the non- Abelian case 
one has to consider the Wilson loop average < W(C) > for the stretched rectangle 
C = r x T with small r and large T: 

V(r) = - lim —In < Tr Pexp (ig I dxaA>iT a ) > . (16) 

T^oo iT V JC J 
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While there is no rigorous way to evaluate < W(C) > analytically, model calculations 
H, [5| do reproduce the same r 2 correction as the leading one at short distances. In 
particular a model of Ref. M gives the following result: 

V(r) = -C F ^- +r 2 $(r) + 0{r 4 ), (17) 
r 

where 

*(r) = ^ < a s (G%) 2 > (V 1 + fr-^Cr)) * (18) 

and p is the characteristic size of non-perturbative fluctuations dominating the gluon 
condensate, < a s {G ( t v ) 2 >. Thus, we have the same r 2 correction in so far as 

<a a [G a f >~ A 4 QCD , p ~ 



3. The examples given above indicate that it is the size of non-perturbative fluctu- 
ations which is crucial to determine the character of leading power corrections to the 
potential at short distances. What unifies renormalons and the model underlying Eq. 
(P~8[), is that in both cases the characteristic size of non-perturbative fluctuations is of 
order Aq Cd . To visualize the connection between the short and large distances in the 
most transparent way let us represent, in an Abelian case, the potential energy of the 
QQ pair as an integral over space from the quark electric fields: 

V(r) = — /dVEi(r') -E 2 (r + r'). (19) 

47T J 

In particular, the Coulomb potential can be obtained of course by integrating over the 
fields E12 of two point charges. On the other hand, if the electric fields are modified at 
large distances, then there arises a correction to the Coulomb energy at small distances 
as well. 

Consider as an example two charges of opposite signs in a cavity of size R. Then 
the electric field of the charges, which is that of a dipole at large distances in empty 
space, changes at r' > R. The corresponding change in V(r) is of order 

which is in agreement with the correction to the static potential discussed above. The 
renormalon- induced correction to the Coulombic potential @ takes into account that 
at r' > Aqq D the fields are distorted by the confinement effects. The corresponding 
correction to the potential is of order 

2 r°° dV „2 * 3 



SV(r) ~ r 2 —g ~ r^ QCD (21) 

QCD V > 

where a s (r') ~ 1 in this estimate. 

From this point of view, it is not at all obvious that the major effect due to con- 
finement is the modification of the Coulomb field of each of the quarks at r' ~ ^qcd- 
Consider a simplified model according to which the electrostatic field of quarks is a 
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correct zeroth-order aproximation only as far as it exceeds some critical value of order 
A 2 

(E a ) 2 ~ A 4 QC D (22) 

while weaker fields do not penetrate the vacumm because of its specific, confining prop- 
erties. From this condition we get an estimate of distances R cr where the electrostatic 
field of quarks is strongly modified: 



< Aqcd (23) 



W 4 

-^6~ - A QCD 



where for simplicity we have neglected the effect of the running of a s (r' 
The change in the potential is then of order: 



2 

limjy ~ ^- ~ a l / 2 rk 2 QCD , (24) 



i.e., we get a linear in r leading correction to the potential at short distances. Note that 
a linear potential in the context of models of a stochastic vacuum has been claimed a 
long time ago |15|], however, that derivation refers to large distances while the estimate 



applies at short distances. It is worth emphasizing, that the estimate fl24|) is 
entirely dependent on the assumption fl22|) as applied to the dipole field of the QQ pair. 
Similarly, the renormalon effect can be visualized as arising from a similar condition 
but applied to a Coulomb-like field of a quark. Since the field of a dipole is weaker 
than that of a charge, the critical value of E 2 is reached at shorter distances and the 
feedback from these distances is stronger. 

It is worth emphasizing that to realize condition (|22| ) with r — > one needs small- 
size nonperturbative fluctuations in the QCD vacuum. Indeed R cr u — ► as r — ► 0. 
This conclusion also fits well with the discussion above. That is why a confirmation of 
the existence of the linear correction (|24D would provide very important insight into the 
mechanism of confinement. If one tries to speculate about what kind of fluctuations 
these could be, it is natural to turn to the dual superconductor picture of confinement 
(for a recent review see [0). Magnetic monopoles are a crucial field configuration 



in this case. The magnetic monopoles of QCD were introduced [18] in the Abelian 



projection of QCD where they appear as singular objects. Although this could be 



an artifact of the gauge fixing [ITH] convincing evidence for existence of monopoles as 



physical objects was obtained in just this gauge (see |L9| and references therein). If the 
physical size of monopoles is indeed vanishing, the linear potential at large distances, 
or the area law for the Wilson loop, could well continue to r — > where it becomes 



a correction to the Coulomb-like potential. Existing data on the lattice [IT] do not 
indicate any change in the linear in r piece of the potential at small distances but no 
special mesurements targeting this behaviour have been performed so far. 

As is mentioned above, very recent measurements do indicate [[Hj the presence of a 



l/Q 2 correction in the effective coupling (see Eq(|])) with a positive constant Cq2\ 

{c Q2 )iat > 0. (25) 
We observe that this sign is in fact opposite to what one would expect from a model in 
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which the effective coupling at large distances is frozen||: 

ttfr * h {lnQVA- 2 Q cD + ^QCD-Q 2 ) ■ (26) 

On the other hand, if one assumes that linear in r potential established at large distances 
continues to short distances as well then the 1/Q 2 correction to the coupling would have 
the same sign as indicated by the data (p5[) . It is worth mentioning once more that 
although most recent discussions [0, |[ emphasize the analyticity aspect of the high-Q 2 
behaviour of the coupling constant the analiticity itself does not fix the form of the 
correction (see, e.g., [pO ). On the other hand, as is discussed above, the existence of 



small-size fluctuations seems to be a prerequisite for the 1/Q 2 corrections. 

4. So far we discussed the static quark potential which can be measured on a 
lattice. Quark potential is also relevant to the physics of bound states of heavy quarks 
(for a review and references see ) . As the zeroth order approximation one considers 



usually Coulomb-like states in the potential (1). Then if one treats the potential K 2 r 
as a perturbation the corrections to the energy levels are: 

SE nl = hsn 2 -l(l + l)]aK 2 , a = — %— (27) 
2 mupa 

where a is the corresponding Bohr radius. Nominally this correction is enhanced by 
powers olm/ Aqcd as compared to correction due to the possible r 2 piece in the potential 
in (see Eq. (HD) or the Voloshin-Leutwyler correction j21 |. In particular, the latter 
correction to energy levels is inversely proportional to m 3 



i 3 - 



, 1.67tt < a s {G a av ) 2 > , 
(SEJvl = TT^ 1 — (28) 



where we have quoted the result for n — 1 which might be the most relevant to phe- 
nomenological applications. There are, however, important caveats to this statement on 
the dominating role of the corrections due to the linear term, if it exists. First, standard 
calculations |22[ already contain a linear term in the potential (extrapolated from large 



distances), so that K 2 shoud be understood rather as a deviation from this form of the 
potential at short distances. Furthermore, the Voloshin-Leutwyler correction contains 
large factors, like a -4 which are important numerically. Since we do not know yet of 
an estimate of the constant K 2 , and also the value of the gluon condensate is not well 
known, it is difficult to compare the relative size of various contributions. The best 
strategy at the moment for the detection of the short- distance K 2 r piece is, as already 
stated, through lattice measurements. 

5. On the theoretical side, the most important reservation about applying and 
comparing various corrections simultaneously is that the potential picture has a limited 
scope of validity. To clarify the applicability of the renormalon-induced and of the 
linear correction (see Eqs. (|9]) and (Q), respectively) to quarkonium physics, let us first 
review the well known results from QED and QCD, 

The particular QED effect which might imitate the interaction of quarks at short 
distances with large scale fluctuations is the shift of atomic levels in a cavity. In the 
static potential picture, the modification of the potential due to the electrostatic inter- 
action of a dipole with the cavity, SV(r) ~ r 2 /R 3 , is found first (see Eq.(^). The energy 
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shifts could be obtained then by averaging this potential over the unperturbed wave 
functions. This procedure would be equivalent to using the r 2 piece in the potential @ 
as a perturbation on the Coulomb-like states. 



A consistent quantum- mechanical treatment of the problem |23[ , on the other hand, 
starts with the dipole interaction: 

H int = -ed-E. (29) 

where E is the electric field associated with zero-point fluctuations of the electromag- 
netic field in the cavity. The spectrum of the zero-point fluctuations depends on the 
size R through the boundary condition, and one is interested in fact only in these in- 
dependent terms. Next, the shift in energy levels E n is obtained to second order in 

Hint'- 

5E n ~ V n k(E„. — Ek + Wchar)" 1 ^. (30) 

We consider bound states and have set V n k = (n\SL in t\k). Since we are interested only 
in estimates, we simply retain the contribution of a characteristic photon frequency, 
^char ~ 1/-R- The result for 5E n depends crucially on the relative magnitude of E n ~ 
ma 2 ~ a/a and c<J c har, i-e., on the relative magnitude of R and a/a, where a is the Bohr 
radius. Note that a new scale, that is a/a, emerges at this stage. Namely, if u chllT ^> E n 
then 



5E ~ aa 2 (E 2 )R ~ a^-r, R<^a/a (31) 

R 6 

since e 2 d 2 ~ aa 2 . The shift corresponds of course to evaluating the matrix element 
of the electrostatic potential, 8V(r) ~ ar 2 / R? discussed above. Thus, in this case the 
potential picture does apply for the evaluation of the energy shifts. 
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On the other hand, if R 3> a/ a then uj c har in the energy denominator of Eq. (|30| ) 
can be neglected and 

3 

5E ~ cra 2 (E) 2 - ~ ^j, R->a/a (32) 
a R* 

where we used E 2 ~ R~ A for the characteristic frequencies. Eq. (^) is in clear violation 
of the potential picture. Eq ( [32]) could be interpreted by saying that the electrostatic 
potential ~ r 2 is replaced, when the distance between the particles a is much smaller 
than aR, by an effective potential ~ r 3 . But this is true only as far as rough estimates 
are concerned. Rigorously speaking, there is no potential whatsoever corresponding to 
the shifts obtained in this way [pi]. Note also that the emergence of the scale R ~ a/ a 



can be understood in a simple way as an effect of retardation. Indeed, the time needed 
to communicate with the distances of order of the size of cavity R can be called the 
retardation time, T ret ~ R. For the potential picture to be valid this time should be 
smaller than the revolution time which is order, T rev ~ a/v ~ a/a. The potential 
picture becomes distorted once T ret ~ T rev . 

In the QCD case one considers pi] atom-like systems of heavy quarks Q with a size 
which is, at least formally, much smaller than the confinement radius, 

a s (M)-M > Aqcd (33) 

where M is the heavy quark mass. In this case the quarks are turning so fast and at 
such a small distance that they cannot be resolved by soft gluons. To account for the 
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interaction of the QQ pair with these soft gluons one starts again with the interaction: 



H mt = - V^(t?-*2)d-E a , (34) 

where, the tf refer to the Q and Q in the quarkonium, and by E a one understands 
the soft gluonic fields (vaccum fields). While the exact form of these fields is an as- 
yet unsolved strong-coupling problem, one assumes usually |] that there is a typical 
frequency, uj c h ar ~ Aq C d- The intensity of the vacuum fields is characterized by a 
vacuum expectation value {(y. s {G a ^ u ) 2 ) which is treated phenomenologically. 

It may be worthwhile here to discuss the situation represented by Eq. ( |3~3"|) from 
the renormalon viewpoint. In order to do this we consider the shift in the interaction 
energy of the quarks in the quarkonium and look for the infrared sensitive contributions 
in perturbation theory. This is seen to arise from the class of diagrams shown in Fig. (la). 
The external states are the color singlet quarkonium states, the exchanged gluon is soft 
(renormalon chain) and the intermediate state in this figure is a continuum color octet 
state which is short lived because of its high virtuality. Indeed the energy denominator 
for this intermediate state is ~ (E — e) ~ B, where E is its energy, and e that of the 
external color singlet quarkonium state with binding energy B . Since B ^> Aq C d, 
(Eq.(|33"D), we see that the intermediate state is far off shell and we thus arrive at 
the reduced diagram shown in Fig (lb). The renormalon chain on the soft gluon line 
in this figure is the standard representation of the ( perturbative ) gluon condensate 
< a s G 2 > which drives this contribution. A standard analysis ( see e.g., E3| ) produces 



a contribution proportional to Aq CD t . It should be emphasized that the situation 
being discussed, i.e., T rev <C Ag^ D ,is precisely the case when the potential picture is 
not valid for the non coulombic corrections. 

Both the cases of static quarks [f|, Q and of atom-like systems (quarkonium) [PTJ] 
with Eq ^> Aqcd were considered in the literature. In the latter case one arrives at a 
natural generalization of Eq. ( |32| ) : 

A 4 

5E ~ra 6 - QCD A m (35) 
{ma s p 

where we have indicated also the sharp dependence on the principal quantum number 

n. Note also that Eq. fl3"5] ) has one extra power of a is the denominator as compared 

2 



with (R!3) because it is the product a s (G a ) which is renormalization group invariant 



It is worth emphasizing also that in the QCD case there exists an extra problem 
with applying the potential picture to the bound QQ states. Namely, the energies of 
intermediate states in Eq. ( |30D are now energies of color states since gluons carry color. 
Formally, these are QQ continuum states. Color states are of course widely used in 
perturbative QCD, and in this respect the situation does not look exceptional. The 
standard constraint, however, is that color states can be excited only for short times 
and Eq. (|30D , when applied in QCD, is therefore consistent as long as the energy 
denominator is much larger than Aqcd, i-e., in the situation discussed above. 

These considerations suggest that, to be rigorous, the only case which can be treated 
consistently is when 

En ~ l " 5 ^ » Aqcd (36) 
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This implies that it is not just the size of the system that has to be smaller than Aqqd, 
but that the binding energy has to be large compared to Aq C d- in the real world, 
where the condition is never satisfied by a wide margin, only a detailed calculation can 
tell the extent and degree of the applicability of the formalism. In this respect it would 
seem that reasonably good results may be obtained for bb with n — 1 and, to a lesser 
extent, n = 2 and cc with n = 1 |T3. 

Returning to the case of the linear correction to the potential (f|) we observe that 
for large enough quark masses, the retardation effects wipe out any kind of a potential 
so that the Voloshin-Leutwyler regime sets in. However, for the linear correction it 
happens at larger masses than for the quadratic correction since it is associated with 
shorter distances R cr ~ cn^a^Ag^ (see fl23|)) where a is the Bohr radius. Therefore, 
the retardation and revolution times get comparable if 

a < o?' A k- Q l CD (37) 

where a is understood to be a function of the mass (and of the quantum number n). 
There is an extra power of a on the right hand side compared to the case of the standard 
correction. 

9. To summarize, we have argued that the assumption that weak color fields do 
not penetrate the vacuum, implies a linear correction to the static quark potential at 
short distances. For this hypothesis to be realized there should exist small-size non- 
perturbative fluctuations. If these fluctuations are monopoles then the linear potential 
at large distances could possibly be extrapolated to short distances which does not con- 
tradict any (lattice) data at the moment. The recent indication [|H| that the running 
coupling has a 1/Q 2 correction does imply a linear correction to the potential as well. 
If confirmed, therefore, it would indicate the existence of small-size non-perturbative 
fluctuations. We have also considered phenomenological implications of the linear cor- 
rection to the potential for quarkonia but found that at this time no conclusion can be 
drawn on the existence or absence of such a correction. 
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Fig. la Class of diagrams contributing to the shift in 
interaction energy of heavy quarks in quarkonium 




Fig. lb Reduced diagram for far off-shell intermediate 
states. 
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